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Abstract 

We give the conditions for a system of N- coupled Korteweg de 
Vries(KdV) type of equations to be integrable. Recursion operators 
of each subclasses are also given. All examples for N = 2 are explicitly 
given. 
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I. Introduction 

In [|l| we gave an extension of the recently proposed Svinolupov Jordan KdV 
[|] systems to a class of integrable multicomponent KdV systems and gave 
their recursion operators. This class is known as the degenerate subclass of 
the KdV system. In this work we will extend it to a more general KdV type of 
system equations containing both the degenerate and non-degenerate cases. 
This is a major step towards the complete classification of KdV systems. In 
addition we give a new extension of such system of equations. 
Let us consider a system of nonlinear equations of the form 

qI = ^qLx + s]kq'ql. (1) 

where i^j,k = 1,2,...,N , are functions depending on the variables x , 
t , and 6* , 5**^ are constants. The purpose of this work is to find the 
conditions on these constants so that the equations in (|l|) are integrable. In 
general the existence of infinitely many conserved quantities is admited as 
the definition of integrability. This implies the existence of infinitely many 
generalised symmetries . In this work we assume the following definition for 
integrability: 

Definition: A system of equation is said to be integrable if it admits a 
recursion operator. 

The recursion operator (if it exists) of the system of equations given in (|^) 
in general, may take a very complicated form. Let the highest powers of the 
operators D and D"^ be respectively defined by m=degree of R and n=order 
non-locality of R. In this work we are interested in a subclass of equations 
admitting a recursion operator with m = 2 and n = 1. Namely it is of the 
form 
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R) = h)D' + a),q'^ + c)^qlD~' (2) 





-I ■ 



where D is the total x— derivative , D is the inverse operator and a*;, , c*;, 
are constants with 

Before starting to the classification of (|I]) we recall a few fundamental prop- 
erties of the recursion operator. An operator i?* is a recursion operator if it 
satisfies the following equation 

R), = Fl^R] - F;^ (4) 
where F^* is the Frechet derivative of the system (|I|) which is given by 

ai = i^'V^ (5) 

where (T*'s are called the symmetries of the system (|I|). The condition (|^) im- 
plies that equation (|l]) is itself is assumed to be in the family of the hierarchy 
of equations (or flows) 



where for all n = 0,1,2,... , o"^ denotes the symmetries of the integrable 
KdV system (|I|). For instance for n = 0, 1 we have respectively the classical 
symmetries cTq = g*^ and cr} = q\. 

Eq.(^ is called the symmetry equation of (P with 



F;^ = b^D' + s),q^ + 4^q'D. (6) 
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Recursion operators are defined as operators mapping symmetries to sym- 
metries , i.e 

Ria^ = \a\ (7) 

where A is an arbitrary constant. Equations and imply (p. It is the 
equation @j which determines the constants 0*^ , c*-;, in terms of 6*- and Sj^'s. 
The same equation (^) brings severe constraints on U- and s*;.. 
We shall obtain a classification of (|I]) based on the matrix 6*-. 

(I) , detih]) =0 , 

(II) . detib]) ^ 0, 

and also we divide the classification procedure, for each class, into two parts 
where s*^ = s\j and s*^ 7^ s\y For the system of equations admitting a 
recursion operator we have the following proposition 

Proposition 1: Let q^{t,x) be functions of t and x satisfying the KdV 
equations ([l|) and admitting a recursion operator Rj in (§) then the constants 
6* , s*;, , a*;, , c*;, satlsfy (in addition to the (^) the following relations 



0, 



0, 
0, 



_L /^^ 



0, 



0, 



'^jk^lm ^km'^jl ^Ik'^jm k'^j m ~^ ^ j I'^km ~^ '^k j m 



0, 



■'kniK 



P3 



3V' 



0. 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 
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Now we will discuss the problem of classifying integrable system of equations 
(|l]) for the two exclusive cases depending upon the matrix by 

II. Classification for the class det{Uj) = 0. 
In this subclass we assume the rank of the matrix 6* as — 1. Investigation 
of the subclasses for other ranks of matrix b can be done similarly. For this 
case we may take 6* = dj — k'^kj where ki is a unit vector , k^ ki = 1. In this 
work we use the Einstein convention , i.e., repeated indices are summed up 
from 1 to . We find the following solution of (PHI^) for the parameters a*- 
and djj. for all 

Proposition 2: Let k^ be a constant unit vector and 6*- = dj — kj then 
the complete solution of the equations are given by 

2.1. 

" 3 + 3 ~ + + b')kiky, 

= ^sij - i [k'{ki ttj + kj ni) + {-ak' + U)kikj\ + k'kiUj, (15) 

where 



Tlj — k kj^S j Qj kjj 

cij —— k kj^Sj^j Ti kjj 

U = k^k^4^^-nk\ (16) 
where a = k"'an, n = k^Ui and s*.;.'s are subject to satisfy the following 



'^jm ^Ik 1 



^jk 



— S 



kj 



F [kj{ak - Uk) - kk{aj - nj)], 



(17) 
(18) 
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kn s'lj = nikj + kittj , (19) 

k^si^ = k% + b'kj, (20) 

r sl^ = {n- a)k'kj + k'ttj + b'kj, (21) 

n [{ai - rii) kj - {aj - rij) h] = 0, (22) 

ai = pni + aki, (23) 



where p is a constant. At this point we will discuss the classification proce- 
dure with respect to the symbol s'jj^ whether it is symmetric or non-symmetric 
with respect to its lower indices. 

(A). The symmetric case, s*;. = s^j : Among the constraints listed in 
proposition 2 the one given in (|2^) implies that s^-^'s are symmetric if and 
only a ai — rii = a ki where a = n — a. There are two distinct cases depending 
on whether n = or n 7^ 0. We shall give these two sub -CcLSGS cLS corollaries 
of the previous proposition. 

Corollary 1: Let 5*^^. = s\j and n = 0. Then we have the following solution 
for all N. 

2 1 

= 34^ + 3 ~ + ^']' 

4 = ^ 4 - ^ Wif^i «fc - 2^fc rij ) + kk kj ¥] , (24) 
where a = , p = 1 and 

ai = ni, ni = hkhij, U = k^k^si^. (25) 
The vector A;* and s*^ are not arbitrary , they satisfy the following constraints 

- sL^m = -2 {kj ni - k rij) (-F + U kj, 
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nikj + ki Tij, 
k^si^ = k% + b'kj,. (26) 

As an illustration we give an example to this case [|l|. A particular solution 
of the equations ( ^61) for = 2 



^5 = - Vi = 
3 

2 

where i,j = 1,2 and 



3 cy. 

^] k = 9 «1 + Ot2x' Vj Vk + -^y' iVj Xk + Vk Xj), (27) 



x' = Sl , y' = 5l (28) 

and 

ki = yi , ni = -aiXi , = a2x\ (29) 
Constants a*;^, and d^^ appearing in the recursion operator are given by 



a] f. = aix' Xj Xk + a2 x' yj yk + ^y' Xj yk, 



c 



■jk 



— X Xj Xk + —y Xj yk- (30) 



2 ' 2 



Taking ai = 2 and a2 = I (without loss of generality) we obtain the following 
coupled system 



vt = {uv)x. (31) 
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The above system was first introduced by Ito 0] and the biHamiltonian 
structure has been studied by Olver and Rosenau . The recursion operator 
of this system is given by 

R= \ \ . 32 

\^ v + v^D-^ 

In we have another example for = 3. 
For the case n for all we have 

Corollary 2: Let 5*^ = s\j , n ^ and p = then the solution given in 
proposition 2 reduces to 



2-1 



3 ^' 3 
1 . 1 



c 



slj --{a-2n)k' kj k (33) 



where 



Uj = nkj, ttj = akj, U = ak\ (34) 
and the constraint equations 



kn s^j = {a + n) kikj 

r s]^ = {a + n) k'kj (35) 



For this case we point out that solution of (^) and ( pS]) gives decoupled 
systems. 
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(B). Non-symmetric case s*;, 7^ s\j: In this case tlie constraints in propo- 
sition 1 , in particular (p2D implies that we must have n = 0. In this case we 
have the following expressions for a^j and c^j for all 

2.1. 

o + o W (kj ai - 2ki Uj) + (-a F + h')kikj\, 

o o 

4, = - I Wih Uj + kj m) + (-a k' + U)kikj\ + k'kmj, (36) 

where 



Qjj — h h-iSj^jj 

n = k^k'sl,. (37) 
and the constraint equations among the parameters are 



+ (^^ - = 0. (39) 

For = 2 we will give an example. Constants a*jr. and djj^ appearing in the 
recursion operator are given by 



^jk 



tti x' Xj Xk + a2 x' Vj Vk + as y' Xj yu, 

«i i , i 

— X XjXk + aiy XjVk- 



(40) 
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where ai,a2 and 03 are arbitrary constants and 

ki = Hi , ni = a^Xi , U = a2x\ = aiXi. (41) 
We obtain the following coupled system 

ut = Uxxx + Saiuux, 

vt = asUxV + aiuvx- (42) 
which is equivalent to the symmetrically coupled KdV system 



Ut = Uxxx + Vxxx + QuUx + AuVx + 2UxV, 

Ut = Uxxx + Vxxx + QvVx + 4:VUx + 2VxU (43) 

and the recursion operator for this integrable system of equations ( ^2]) is 

/ D^ + ai(2u + UxD-^) \ 
R= \ ^ _ . (44) 

III. Classification for the case det{b^j) ^ 0. 
As in the degenerate case det{b^j) = we have two sub-cases, symmetric and 
non- symmetric, before these we have the following proposition 
Proposition 3: Let det{b'j) ^ then the solution of equations given propo- 
sition 1 is given as follows 



^i, = \i3^i,-^s),)^\c\h^^Kf^. (45) 
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where 



Kh = ^,-s',i (46) 

and the constraint equations 

b'ic],-U,c], = Q, (47) 

bCTK^ - ClK"^ - C^Kl^l = (48) 

c'^^Ki, + = (49) 

^)k^hn ~ ^]m^\k = (^O) 

Kt<^ = 0, (51) 



i'^]k ^)sf„ + (c^„ al^)sfj + Kl^f^a^ji + K^^a\i — 0. (52) 
where is the inverse of 6[. 

(A). Non-symmetric case s*;, ^ s\j : Eqs. p7| - |52| ) define an over- 
determined system for the components of s^j^, . Any solution of this system 
leads to the determination of the parameters a*^ and c*; by (^51) . 
As an example we give the following,for N = 2, coupled system 

Ut = -iaUxxx + ^bUxV + buVx, (53) 

where a and b are arbitrary constants. This system ,under a change of vari- 
ables, is equivalent to the KdV equation with the time evolution part of its 
Lax equation [0. The recursion operator of the system (|53| ) is 
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R 



^ |(3aD2 + bv) |(3m + 2u^D-^] 



(54) 



\^ ^{3aD^ + Abv + 2bv^D- 

Hence the KdV equation coupled to time evolution part of its Lax-pair is 
integrable and its recursion operator is given above. This is the only new 
example for = 2 system. 

(B). Symmetric case s^j^. = s^^- : When the symbol s*;^. is symmetric 
with respect to sub-indices the parameters Kjj^ vanish. Then the equations 



(^5| - p^ ) reduce to 



2 ■ 1 • 

where the parameters b]. and s*^ satisfy 

&Hfc-&14 = 0, (56) 

^]k^tm ~ ^Ik^^m — ^- (57) 

We shall not study this class in detail , because in [|l|| some examples of this 
class are given for N = 2. Here we give another example which correspond 
to the perturbation expansion of the KdV equation. Let = 5* u where 
i = 0,1, 2, N and u satisfies the KdV equation u^t = u^^xx + Guu^^. The 
g*'s satisfy a system of KdV equations which belong to this class 



c^!t = q%x + (^qyx, (58) 

9i = 9^ + 6 (59) 
^ = 9f... + 6[(?T + gN'],x, (60) 
(61) 



12 



9? = 9j.. + 3E[^^(ffl,^ (62) 

i=l 

IV. Fokas-Liu Extension 

The classification of tlie KdV system given in tliis work witli respect to tlie 
symmetries can be easily extended to the following simple modification of 

m 

4 = ^qL. + s'jkQ'Q^. + x) q%, (63) 

where Xj's are arbitrary constants. Eg. (p3|) without the last term will be 
called the principle part of that equation. Hence the equation (|lD we have 
studied so far is the principle part of its modification (|63D. We assume the 
existence of a recursion operator corresponding to the above system in the 
form 



R^ = b^D' + a^,q' + c], ql D'' + w] (64) 

where w^jS are constants. We have the following proposition corresponding 
to the integrability of the above system 



Proposition 4: The operator given in (|6^) is the recursion operator of the 
KdV system ( |63D if in addition to the equations listed in proposition 1 
the following constraints on the constants and are satisfied 



Xi a]k - X\ 4 - X\ Cji + Xj 4/ - sli + < ' 



0, 



Xf '^ji ~ 0' 

(x? - ^ • ) 4 - (xl -K)b^ = 0, 



(65) 
(66) 
(67) 
(68) 
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Xj all - Xl a% + < sfj - w] s],, = 0. (69) 

Since the constraints (Q-jl^) are enough to determine the coefficients a*;, and 
dji^ with some constraints on the given constants 6* and s*;,'s we have the 
following corollary of the above proposition 

Corollary 1: The KdV system in ( |63D is integrable if and only if its principle 
part is integrable. 

The principle part of (^) is obtained by ignoring the last term (the term with 
Xj)- Hence the proof of this corollary follows directly by observing that the 
constraints on the constants x] and Wj listed in (|65|-p9|) are independent of 
the constraints on the constants of the principle part listed in (^0). Before 
the application of this corollary let us go back to the proposition 4 and ask 
the question whether the KdV system (|63D admits a recursion operator with 
w'j = 0. 

Corollary 2: The KdV system (^3]) admits a recursion operator of the 
principle part then the last term x] ^ symmetry of the principle part. 
If Wj = the above equations ( |65| - |69|) reduce to 



Xi a% - xi 4 - Xi 4 + Xj 4z = 0, (70) 

X^s],-xXs^, = 0, (71) 

Xf bl - Xfc = 0, (72) 

Xf c], - x", 4 = (73) 

In order that the term x] be a symmetry of the principle part the 

constants x] are subject to satisfy the following equations 



Xf s)k - X\ 4 = 0' 



(74) 
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X-ftfc-Xfcfe) = 0, (75) 
xfK;, + x^Ki, = 0, (76) 

These equations simply follow from the set of equations ( |70| - |73D and hence 
the quantities a* = Xj Q.^x symmetries of the principle part. 

By the application of corollary 1, the full classification of the system (|6^) 
with the recursion operator ( |64|) such that q\ (i.e., the system of equations 
themselves) belong to the symmetries of this system is possible. To each 
subclass given in the previous sections there exists a Fokas-Liu extension 
such that * = x] with the following constraints 

Xl4-xfc}, = 0, (77) 
Xz' 4fc - 4 < - (4 - = 0. (78) 

The above constraints are identically satisfied for the class {det{b'j) ^ 
, symmetrical case) when x*- = a5* + Hence the Fokas-Liu extension 
of non-degenerate-symmetrical case is straightforward with this choice of x*. 
Here a and /3 are arbitrary constants. 

For the degenerate case the set of equations (|77| - |78|) must solved for a 
given principle part , fe* and a*;.. Recently a system of integrable KdV system 
with N = 2 has been introduced by Fokas and Liu . This system is a nice 
example for the application of corollary 1. We shall give this system in its 
original form first and then simplify 



u,t + v^,, + (3 A + 2 {3^) u M,, + (2 + /3i (3^) (uv),, + 

Pi Ps V v^^ + (/5i + (3i) (32 ^i.xxx + (1 + A Pa) P2 v,xxx = 0, (79) 
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v^t + + (2 + 3 A [3i) V v,^ + {pi + 2 /94) (uv),^ + 
/3i /?3 /34 + + Pa) P2 Pa m,..^ + (1 + A /?4) /34^^,xxx = 0. (80) 

where Pi, P2 , Ps , and P4 are arbitrary constants. The recursion operator of 
this system is given in p|. Consider now a hnear transformation 

u = mi r + riis , v = 1712 r + 77,2 s (81) 

where ttii , 1122 , ni , 72,2 are constants and s , r are new dynamical variables, 
= (s,r). Choosing these constants properly, the Fokas-Liu system reduces 
to a more simpler form 



r,t = (rs),a: + air^ + 02 Sx, 
s,t = 7i s,xxx + ^2rr^^ + 3s s^x + "3 + "4 s,^, (82) 

where we are not giving the coefficients ai, 02,013 and 04 in terms of the 
parameters of the original equation given above, because these expressions 
are quite lengthy. The only condition on the parameters is given by 
«3 = 72 «2- this guarantees the integrability of the above system (|82D . On 
the other hand the transformation parameters are given by 



I3i + Pa . 
m2 = - mi , n2 = PAni, (83) 

i + Pl P4 

ni = -T^ , S = Pi{l + Pl) + 2p,. (84) 



The principle part of the Fokas-Liu system ( ^2]) is exactly the Ito system 
given in (pl|) and hence the recursion operator is the sum of the one given in 
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2|) and x) which are given by x\ = , X2 = «3 , X? = «2 , X2 = "i- 



That is 



R 



(85) 



-fi D"^ + 2 s + s^D ^ + 04 72r + a3 
Another example is given very recently in very different context for N = 2 
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u,t = - v^^x^ + 2uv^^ + u^^ V, 
V t = Svv X + 2 au X- 



(86) 



The principle part of this equations is transformable to the Fuchssteiner 



system given in (|4^) . Taking ai = 2, 03 = 4 and scaling x and t properly we 
get (without loosing any generality) 



i"t 2 Txxx ~\~ Grr^;, 

St = 2sxr + Asrx. (87) 

The transformation between the principle part (a = 0) of ( PUf ) and ( P7D is 
simply given hj u = nii s + |r,f = 2r. Then the recursion operator of the 
system (^) is given by 

y2a 2v + VxD-^ j 

V. Conclusion 

We have given a classification of a system of KdV equations with respect to 
the existence of a recursion operator. This is indeed a partial classification. 
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Although we have found all conditions for each subclassess we have not pre- 
sented them explicitly. We obtained three distinct subclasses for all values of 

and gave the corresponding recursion operators. We also gave an exten- 
sion of such systems by adding a linear term containing the first derivative of 
dynamical variables. We called such systems as the the Fokas-Liu extensions. 
We proved that these extended systems of KdV equations are also integrable 
if and only if their principle parts are integrable. For A = 2 , we have given 
all subclassess explicitly. Among these the recursion operator of the KdV 
coupled to the time evolution part of its Lax-pair seems to be new. Here 
we would like to add that when N = 2 recursion operators , including the 
Fokas-Liu extensions, are hereditary ||1CI|| . 

Our classification crucially depends on the form of the recursion operator. 
The recursion operators used in this work were assumed to have degree two 
(highest degree of the operator D in R) and non-locality order one (highest 
degree of the operator in R). The next works in this program should be 
the study on the classification problems with respect to the recursion opera- 
tors with higher degree and higher non-localities. For instance when N = 2 , 
Hirota-Satsuma, Boussinesq and Bogoyavlenskii coupled KdV equations ad- 
mit recursion operators with m = 4 and n = 1 |1T[. Hence these equations 
do not belong to our classification given in this work. 

This work is partially supported by the Scientific and Technical Research 
Council of Turkey (TUBITAK). M.G is a member Turkish Academy of Sci- 
ences (TUBA). 
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